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Abstract 

A measurement is presented of dijet and 3-jet cross sections in low-|f| diffractive deep- 
inelastic scattering interactions of the type ep — > eXY, where the system X is separated by 
a large rapidity gap from a low-mass baryonic system Y. Data taken with the HI detector 
at HERA, corresponding to an integrated luminosity of 18.0 pb _1 , are used to measure 
hadron level single and double differential cross sections for 4 < Q 2 < 80 GeV 2 , xp < 
0.05 and prjet > 4 GeV. The energy flow not attributed to jets is also investigated. 
The measurements are consistent with a factorising diffractive exchange with trajectory 
intercept close to 1.2 and tightly constrain the dominating diffractive gluon distribution. 
Viewed in terms of the diffractive scattering of partonic fluctuations of the photon, the data 
require the dominance of qqg over qq states. Soft colour neutralisation models in their 
present form cannot simultaneously reproduce the shapes and the normalisations of the 
differential cross sections. Models based on 2-gluon exchange are able to reproduce the 
shapes of the cross sections at low xjp values. 
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1 Introduction 



The observation of deep-inelastic scattering (DIS) events at HERA containing a large gap in 
the rapidity distribution of the final state hadrons [[]]] has generated considerable renewed in- 
terest in understanding colour singlet exchange in strong interactions. At high energy, such 
interactions are interpreted as being due to diffractive scattering. HERA has made it possible to 
study diffraction using a highly virtual photon probe. This offers the chance to illuminate the 
underlying dynamics in terms of quantum chromodynamics (QCD). 

Inclusive diffractive DIS is principally sensitive to the role of quarks in the scattering pro- 
cess More insight into the gluonic degrees of freedom can be obtained by studying the 
hadronic final state Final states containing heavy quarks or high transverse momentum 
(pr) jets are of particular interest, since the additional hard scales may ensure the applicability 
of perturbative QCD techniques p-|12|]. High p T jet production in diffraction has previously 



been studied both in pp collisions [|T3|-|T6]] and at HERA [|7],§p. 

In this article, a high statistics measurement of diffractive jet production is presented, which 
was performed using the HI detector. The data were obtained using events where the proton (or 
a low-mass proton excitation) loses only a small fraction of its incoming momentum and escapes 
undetected through the beam pipe. Separated from this system by a large rapidity region devoid 
of activity, the hadronic system X is well contained within the central part of the detector and 
contains the high p T jets. The luminosity is increased by an order of magnitude compared with 
previous HI measurements ^ and the kinematic range is also extended. This makes it possible 
to extract double differential cross sections for the first time and to study 3-jet as well as dijet 
production. 

The dijet data yield direct constraints on the diffractive gluon distribution and are used 
to investigate the QCD [JT7p and Regge [ |T%| ] factorisation properties of diffractive DIS. QCD 
inspired models [[IT], [L2|, |19|] based on the exchange of a pair of gluons from the proton [ |20| ] 
are compared with the data in a restricted kinematic region where they are most likely to be 
applicable. Predictions from soft colour neutralisation models pT[-P3|] are also confronted with 
the data. 

The article is organised as follows. The kinematics of diffractive scattering at HERA are 
introduced in section In section [3J, an overview of phenomenological models and QCD cal- 
culations relevant for diffractive jet production is given and the Monte Carlo simulation of 
diffractive events is described. In section ^, the data selection, the cross section measurement 
procedure and the determination of the systematic uncertainties are explained. The results, ex- 
pressed in terms of hadron level single and double differential cross sections, are presented and 
discussed in section |5|. 



2 Diffractive Scattering at HERA 



i 




Figure 1: The generic diffractive process at HERA, where the electron (A;) couples to a photon 
(q) which interacts with the proton (P) via net colour singlet exchange, producing two distinct 
final state hadronic systems X and Y. If the masses of X and Y are small compared with W, 
the two systems are separated by a large gap in rapidity. 

2.1 Inclusive Diffractive Scattering 



Fig. [1] illustrates the generic diffractive process at HERA of the type ep — > eXY. The electron|] 
(with 4-momentum k) couples to a virtual photon 7* (q) which interacts with the proton (P). 
The usual DIS kinematic variables are defined as 



Q 2 



-Q 



y 



P-q 
P-k 



x 



-r 



2P-q 



(1) 



The squared invariant masses of the electron-proton and photon-proton systems s and W 2 are 
given by 



s = (k + Pf ~ Q 2 /xy ~ (300 GeV) : 



(g + P) 2 ~y S -Q 2 



(2) 



If the interaction takes place via colour singlet exchange, the photon and proton dissociate to 
produce distinct hadronic systems X and Y, with invariant masses M x and M Y respectively. 
In the case where M x and M Y are small compared with W, the two systems are separated by 
a large rapidity gap. The longitudinal momentum fraction xp of the colourless exchange with 
respect to the incoming proton and the squared four-momentum transferred at the proton vertex 
t are then defined by 

q-(P- p Y ) 



Xp 



q-P 



t = (P 



Vvf, 



(3) 



where p Y is the 4-momentum of Y. In the analysis presented here, t and M Y are not measured 
and hence are integrated over implicitly!]. In addition, the quantity (3 is defined as 



Xjp 



Q 2 



2q-(P-p Y ] 



(4) 



In an interpretation in which partonic structure is ascribed to the colourless exchange, (3 is the 
longitudinal momentum fraction of the exchange that is carried by the struck quark, in analogy 
to x in the case of inclusive scattering. 

'From now on, the word 'electron' will be used as a generic term for electrons and positrons. 
2 It is noted that for this analysis My — M p dominantly. 
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Proton rest frame 



(a) 



Y* 




(c) 



Y* 



222 




Proton infinite momentum frame 

(b) 




Figure 2: Diffractive scattering in the proton rest frame and the proton infinite momentum 
frame (figure after [pT]|). In the proton rest frame, the virtual photon dissociates into a qq 
state (a), scattering off the proton by colour singlet (e.g. 2-gluon) exchange. In the infinite 
momentum frame, this can be related to diffractive quark scattering (b). The emission of an 
additional gluon forms an incoming qqg state (c). If the gluon is the lowest pt parton, this 
contribution can be related to diffractive Boson-Gluon-Fusion (d). 

2.2 Diffractive Dijet Production 

Viewing DIS at low x in the proton rest frame, the virtual photon splits into a qq pair well 
in advance of the proton target (fig. ga). The qq state may then scatter elastically with the 
proton. The production of high p T final states by the diffractive qq scattering process is heavily 
suppressed [^5p and the invariant masses Mx produced are typically small. It is thus expected 
that for large values of Mx or pr, 0(a s ) contributions due to the radiation of an extra gluon 
become important [f5[^5[]. The result is an incoming qqg system (fig. ||c). 



In the proton infinite momentum frame, the lowest order (i.e. 0(a Q s )) contribution to the 
diffractive cross section is the quark scattering diagram (fig. 0b). The 0(a s ) contributions are 
Boson-Gluon-Fusion (BGF) and QCD-Compton (QCDC) scattering. Unlike inclusive diffrac- 
tive scattering, jet production is directly sensitive to the role of gluons in diffraction due to the 
direct coupling to the gluon in the BGF diagram (fig. ||d). 

There is a correspondence between the proton rest frame and the infinite momentum frame 
pictures, which is discussed here in the context of the leading log(Q 2 ) approximation. For the 
dominant configuration in which the photon longitudinal momentum is shared asymmetrically 
between the partons, diffractive qq scattering (fig. |2|a) can be related to the diffractive quark 
scattering diagram (fig. ^b). If the gluon is the lowest pt parton, the diffractive scattering of 
asymmetric qqg configurations (fig. §c) can be related to diffractive BGF (fig. [|d). If the q or q 
is the lowest p T parton, the process corresponds to diffractive QCDC scattering (not shown). 



Using the non-zero invariant mass squared s of the two highest p T partons emerging from 
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the hard interaction in the 0(a s ) case, the quantity zjp is introduced: 

z P = P-(l + s/Q 2 ) . (5) 

Similarly to j3 for the case of the lowest order diagram (fig. [2]b), zjp corresponds to the longitu- 
dinal momentum fraction of the exchange which takes part in the hard interaction (fig. ||d). 



3 Phenomenological Models and Monte Carlo Simulation 

In this section, several phenomenological approaches and QCD calculations are discussed, 
which attempt to describe diffractive DIS, including diffractive jet production. The focus is 
on the models which are compared with the data in section |[ 



3.1 Diffractive Parton Distributions 

In the leading log(Q 2 ) approximation, the cross section for the diffractive process 7*p — > p'X 
can be written in terms of convolutions of universal partonic cross sections a 1 * 1 with diffractive 
parton distributions ff, representing probability distributions for a parton i in the proton under 
the constraint that the proton remains intact with particular values of xjp and t. Thus, at leading 
twist, 



dxp dt = ^ Jx ^' ^ * ^ ' XJP ' ' ' 



(6) 



This factorisation formula holds for large enough Q 2 and fixed x, xjp and t. This ansatz, in- 
troduced in ^7p, was rigorously proven for inclusive diffractive lepton-hadron scattering 
in fli7|,[28l]. The diffractive parton distributions are not known from first principles, though they 
should obey the DGLAP [E9p evolution equations. Recently, there have been attempts to cal- 



culate the diffractive parton distributions at a starting scale /ijj for QCD evolution under certain 
assumptions. In J30|] , the proton is replaced by a small-size pair of heavy quarks, such that per- 
turbation theory can be applied. A different approach is the semiclassical model by Buchmuller, 
Gehrmann and Hebecker [^Ij], based on the opposite extreme of a very large hadron. In spite 
of the different assumptions, the two approaches give rather similar results for the diffractive 
parton distributions. The general behaviour is the same as the momentum fractions tend to or 
1 and the gluon distribution dominates. 



3.2 Resolved Pomeron Model and Pomeron Parton Distributions 

The application of Regge phenomenology of soft hadronic high energy interactions to the con- 
cept of diffractive parton distributions (section |3T| ) leads to the Ingelman-Schlein model of a 
'resolved pomeron' with a partonic structure [ |Tg| ] invariant under changes in xp and t. The 
diffractive parton distributions then factorise into a flux factor f F / p and pomeron parton distri- 
butions ff: 

f°{x, Q 2 , xjp, t) = fjp /p (xjp, t) ■ ff(/3 = x/xjp, Q 2 ) . (7) 
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The universal flux factor describes the probability of finding a pomeron in the proton as a 
function of xjp and t. The pomeron parton distributions are usually expressed in terms of (3. 

The triple differential cross section for inclusive diffraction d 3 a/d(3 dQ 2 dxp is often pre- 
sented in the form of a diffractive structure function F^ 3 \(3, Q 2 , xp). In the HI collabo- 
ration interpreted a measurement of in terms of a resolved pomeron model: At the largest 
xp studied, it was necessary to consider more generally contributions from sub-leading reggeon 
exchanges! as we ^ as me pomeron, such that (neglecting possible interference terms) 

f t D (x,Q 2 ,x P ,t) = /p/ P (x P ,i)-/f(^Q 2 ) + /fi/ P (^i)7f(AQ 2 ). (8) 

The flux factors for the pomeron and reggeon exchanges were parameterised in a Regge-inspired 
form: 

f{F,K}/p(x F ,t) = C { p M a;jp l-2- {J P. JR} (*) e 6 { r,K}t ; (9) 

with o;{2P,jR}(t) = «{jp,k}(0) + a'r F M yt. From fits in which the parton densities evolve accord- 
ing to the DGLAP equations, parameterisations of the pomeron quark and gluon distributions 
and values for the trajectory intercepts ap(0) and ap(0) were obtained. The resulting value of 
ajp(O) = 1.203 ± 0.020 (stat.) ± 0.013 (syst.) ± 0.030 (model) is significantly higher than 
that obtained from soft hadronic interactions, where ap(0) — 1.08 pT|]. The parton densities 
extracted for the pomeron are dominated by gluons, which carry 80 — 90% of the exchanged 
momentum throughout the measured Q 2 range. 



3.3 Colour Dipole and 2-Gluon Exchange Models 

In the proton rest frame, diffractive DIS is often treated by considering the qq and qqg photon 
fluctuations (fig. |2|a,c) as (effective) colour dipoles. The diffractive YP cross section can be 
factorised into a squared effective photon dipole wave function and a squared 'dipole cross 



section' for the scattering of these dipoles off the proton [p2|,|33|]. The gross features of the 
diffractive (3 distribution can be deduced from a knowledge of the partonic wave functions of 
the photon alone. According to a recent QCD motivated parameterisation [J34p, longitudinally 
and transversely polarised qq states dominate at high and medium values of (3 respectively, 
whereas the qqg state originating from transversely polarised photons is dominant at low (3. 

Investigating diffractive final states with varying p T probes the dipole cross section as a 
function of the dipole size. Large size, low p T configurations interact with the proton simi- 
larly to soft hadron-hadron scattering. Small size, high dipole configurations lead to hard 
scales which encourage a perturbative QCD treatment of the dipole cross section. The precise 
dynamics of the dipole cross section are not known a priori. However, the simplest realisation 
of a net colour singlet exchange at the parton level is a pair of gluons with cancelling colour 
charges [p0[]. We focus below on two recent colour dipole models [ |TT1 , ^2| , p^ ] based on 2-gluon 



exchange, where the cross section is related to the square of the unintegrated gluon distribution 
of the proton T{x, k\) []35|]. Here, k T is the parton transverse momentum relative to the proton 
direction. Other colour dipole approaches can be found in [|9|,p^|,P6]]. 



throughout this paper, the term 'reggeon' (M) will be used to describe this contribution. 
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The dipole approach has been employed in the 'saturation' model by Golec-Biernat and 
Wiisthoff |j9j|. Here, an ansatz for the dipole cross section is made which interpolates between 
the perturbative and non-perturbative regions of cr 7 * p . This model is able to give a reasonable 
description of F 2 (x,Q 2 ) at low x, which determines the three free parameters of the model. 
The parameterised dipole cross section can be re-expressed in terms of F{x, kj), such that 
the diffractive cross section is predicted at t = 0. Introducing an additional free parameter 
B = 6.0 GeV~ 2 to describe the t dependence as e Bt , the diffractive structure function F 2 °^ is 
successfully described. The calculation of the qqg cross section is made under the assumption 
of strong k T ordering of the final state partons (leading log(Q 2 ) approximation), corresponding 

to4 s) <4 9>9) . 

Cross sections for diffractive qq and qqg production by 2-gluon exchange have been cal- 
culated by Bartels, Ewerz, Lotter and Wiisthoff (qq) flTTl ] and by Bartels, Jung, Kyrieleis and 
Wiisthoff (qqg) [|T2|]. The derivative of the next-to-leading order (NLO) GRV gluon parame- 
terisation [ p7| ] is used for J-(x, k^). The calculation of the qqg final state is performed in the 
leading log(l//3), leading log(l/xp) approximation, such that configurations without strong 
k T ordering are included. The calculations require all outgoing partons to have high p T and are 
thus not suited to describe F 2 D ^ . The minimum value p™* for the final state gluon transverse 
momentum is a free parameter which can be used to tune the model to the overall dijet cross 
section. As for the saturation model, the calculation yields predictions at t = 0. The extension 



to finite t is performed using the Donnachie-Landshoff elastic proton form factor [p8[]. The sum 
of the qq and qqg contributions in this model is hereafter referred to as 'BJLW'. 



3.4 Soft Colour Neutralisation Models 

An alternative approach to diffractive DIS is given by soft colour neutralisation models, which 
naturally lead to very similar properties of inclusive and diffractive DIS final states. In the Soft 



Colour Interaction (SCI) model by Edin, Ingelman and Rathsman Q22|], the hard interaction in 
diffractive DIS is treated identically to that in inclusive DIS. Diffraction occurs through soft 
colour rearrangements between the outgoing partons, leaving their momentum configuration 
unchanged. If two colour singlet systems are produced by such a mechanism, the hadronic 
final state can exhibit a large rapidity gap. In the original SCI model, diffractive final states 
are produced using only one free parameter, the universal colour rearrangement probability, 
which is fixed by a fit to F 2 Dt * z \ The model has been refined recently [[23]] by making the colour 



rearrangement probability proportional to the normalised difference in the generalised areas of 
the string configurations before and after the rearrangement. 



The semiclassical model, which was already mentioned in section pTT], is a non-perturbative 
model based on the dipole approach. Viewed in the proton rest frame, qq and qqg photon 
fluctuations scatter off a superposition of soft colour fields associated with the proton. Those 
configurations which emerge in a net colour singlet state contribute to the diffractive cross 
section [|2~5|]. Assuming a specific model for the proton wave functional pl|], the results are 



formulated as a parameterisation of t-integrated diffractive parton distributions p9|], which are 
determined from a combined four parameter fit to F 2 and F 2 D ^ at low x and xp. 
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3.5 Monte Carlo Simulation 



Monte Carlo simulations are used to determine the corrections to be applied to the data to com- 
pensate for the limited efficiencies, acceptances and resolutions of the detector. The generated 
Monte Carlo events are passed through a detailed simulation of the HI detector and are sub- 
jected to the same reconstruction and analysis chain as the data. 



The main Monte Carlo generator used to correct the data is RAPGAP 2.08/06 pOQ. Events 
are generated according to a resolved (partonic) pomeron model (section Contributions 
from pomeron and reggeon exchanges are included neglecting any possible interference effects. 
The parameterisations of the pomeron and reggeon flux factors and parton distributions are 
taken from the HI analysis of F®^ 1 [BJ]. The pomeron and reggeon trajectories and slope 
parameters (eq. g) are a P (t) = 1.20 + 0.26t, bp = 4.6 GeV" 2 , a R (t) = 0.50 + 0.90t and 
bp = 2.0 GeV -2 . The pomeron parton distributions are taken from the 'flat gluon' (or 'fit 2' 
of [Bp) solution in the leading order DGLAP fits to F^ 3 \ Those of the meson are taken from 



fits to pion data [£?]]]. The renormalisation and factorisation scales are set to /i 2 = Q 2 + p\, 
where pr is the transverse momentum of the partons emerging from the hard scattering relative 
to the collision axis in the ^*p centre-of-mass framef]. The parton distributions are convoluted 
with hard scattering matrix elements to 0(a s ). Intrinsic transverse momentum of the partons 
in the pomeron p2| ] is not included. Charm quarks are produced in the massive scheme via 
Boson-Gluon-Fusion. For the production of light quarks, a lower cut-off in p\ is introduced 
in the 0(a s ) QCD matrix elements to avoid divergences. Higher order QCD diagrams are 
approximated with parton showers in the leading log(Q 2 ) approximation (MEPS) [ |43| ] or the 
colour dipole approach^] (CDM) as implemented in ARIADNE [|4"5|]. Hadronisation is 



simulated using the Lund string model in JETSET []4~6|]. QED radiative effects are taken into 



account via an interface to the HERACLES program p7|] . 



The RAPGAP simulation includes a contribution of events where the virtual photon 7* is 
assigned an internal partonic structure. The resolved virtual photon is parameterised according 
to the SaS-2D set of photon parton densities, which has been found to give a reasonable 
description of inclusive dijet production at low Q 2 



Monte Carlo generators are also used to compare the measured hadron level cross sections 
with the predictions of the phenomenological models and QCD calculations described in sec- 
tions PT2"| - p.4| . All of the predictions are made to leading order of QCD. Unless otherwise stated, 
higher order QCD effects are approximated by initial and final state parton showers. RAPGAP 
is used to obtain the predictions of the resolved pomeron model with different pomeron intercept 
values and parton distributions. It also contains implementations of the saturation, semiclassi- 
cal and BJLW models. Both versions of the SCI model are implemented in the LEPTO 6.5.2/3 
generator [ |50| ] . 



4 This frame is also called the 'hadronic centre-of-mass frame'. 

5 'Colour dipole approach' as an approximation to higher order QCD effects should not be confused with the 



'Colour dipole models' introduced in section 3.3 
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4 Experimental Procedure 



The analysis presented in this article is based on HI data taken in the years 1996 and 1997, when 
HERA collided E e = 27.5 GeV positrons with protons of E p = 820 GeV. The data correspond 
to an integrated luminosity of 18.0 pb _1 . A detailed description of the measurement can be 
found in Q5T|] . This section begins with a short overview of the HI detector, after which the 
data selection is described. Then, the cross section measurement and the determination of the 
systematic errors are explained. 

4.1 HI Detector 

The HI detector is described in detail elsewhere [p2[]. Here, we give a brief description of the 
detector components most relevant for the present analysis. The z axis of the HI coordinate 
system corresponds to the beam axis such that positive z values refer to the direction of the 
outgoing proton beam, often called the 'forward' direction^. 

The interaction region is surrounded by the tracking system. Two large concentric drift 
chambers (CJC), located within a solenoidal magnetic field of 1.15 T, measure the trajectories 
of charged particles and hence their momenta in the range —1.5 < t] < 1.5. The resolution is 
<j{j>t)Ivt — O-Olpr/GeV. Energies of final state particles are measured in a highly segmented 
Liquid Argon (LAr) calorimeter covering the range —1.5 < r] < 3.4, surrounding the tracking 

detectors. The energy resolution is a{E)/E ~ 11%/y^/GeV for electromagnetic showers 

and <j{E)/E ~ 50%/ \J E / GeV for hadrons, as obtained from test beam measurements. The 
overall hadronic energy scale of the LAr is known to 4%. The backward direction (—4.0 < 
r\ < —1.4) is covered by a lead / scintillating fibre calorimeter (SPACAL) with electromagnetic 
and hadronic sections. The energy resolution for electrons is a(E)/E ~ 10%/ JeJg^V. The 
energy scale uncertainty is 0.3% for electrons with E' e = 27.5 GeV and 2.0% at E' e = 8 GeV. 
The electron polar angle is measured to 1 mrad. The energy scale of the SPACAL is known to 
7% for hadrons. In front of the SPACAL, the Backward Drift Chamber (BDC) provides track 
segments of charged particles with a resolution of a(r) = 0.4 mm and ra{<p) = 0.8 mm. The 
ep luminosity is determined with a precision of 2% by comparing the measured event rate in a 
photon tagger calorimeter close to the beam pipe at z = —103 m with the QED Bremsstrahlung 
(ep — > ep7) cross section. 

To enhance the sensitivity to hadronic activity in the region of the outgoing proton, the 
Forward Muon Detector (FMD) and the Proton Remnant Tagger (PRT) are used. The FMD is 
located at z = 6.5 m and covers the pseudorapidity range 1.9 < r\ < 3.7 directly. Particles 
produced at larger r\ can also be detected because of secondary scattering with the beam-pipe. 
The PRT, a set of scintillators surrounding the beam pipe at z = 26 m, can tag hadrons in the 
region 6.0 < 7] < 7.5. 

6 This direction corresponds to positive values of the pseudorapidity r\ = — In tan 9/2. 
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4.2 Data Selection 



DIS events are triggered by an electromagnetic energy cluster in the SPACAL in coincidence 
with a CJC track. Scattered electron candidates are then selected with E' e > 8 GeV in the 
angular range 156° < 9' e < 176°. Various cuts are applied on these candidates in order to select 
electrons and reject background from photons and hadrons. Among these are requirements on 
the width of the shower, the containment within the electromagnetic part of the SPACAL and 
the existence of an associated track segment in the BDC. The z coordinate of the reconstructed 
vertex is required to lie within ±35 cm (± ~3cr) of the nominal interaction point. To suppress 
events with initial state QED radiation, the summed E — p z of all reconstructed final state 
particles including the electron^ has to be greater than 35 GeV. The DIS kinematic variables 
are calculated from the polar angle and energy measurements of the scattered electron: 

Q 2 = 4£ e £>s 2 |; y = 1 - § sin 2 | . (10) 

Z £j e Z 

Events which fulfil 4 < Q 2 < 80 GeV 2 and 0.1 < y < 0.7 are selected. 

The selection of diffractive events is based on the absence of hadronic activity in the out- 
going proton region. No signal above noise levels is allowed in the FMD and PRT detectors. 
The most forward part (jj > 3.2) of the LAr calorimeter has to be devoid of hadronic clusters 
with energies E > 400 MeV. This selection ensures that the photon dissociation system X is 
well contained within the central part of the HI detector and is separated by a large rapidity 
gap covering at least 3.2 < 77 < 7.5 from the Y system. The upper limit in i] implies that the 
Y system escapes undetected through the beam pipe and imposes the approximate constraint 
M Y < 1.6 GeV and \t\ < 1.0 GeV 2 . 

The X system is measured in the LAr and SPACAL calorimeters together with the CJC. It 
is reconstructed using a method that combines calorimeter clusters and tracks whilst avoiding 



double counting [p3Q. The dissociation mass is then calculated according to 

Mi = (E,^) 2 -(E,P J ) 2 , (ID 

where the sum runs over all reconstructed objects except for the scattered electron^. W 2 is 
calculated according to eq. (Q). x F and (3 are then computed from 

_ Q 2 + M 2 X Q _ Q 2 

Xjp -g^Tl^ ; q^~+mJ ' (12) 

A cut xp < 0.05 is applied to suppress contributions from non-diffractive scattering and sec- 
ondary exchanges. The resolution in logxjp is approximately 8%. 

The 4-vectors of the hadronic final state particles associated to the X system are Lorentz- 
transformed to the YP centre-of-mass frame, where they are subjected to the CDF cone jet 



algorithm [ |54| ] with a cone radius of J (Ar/) 2 + (A0) 2 = 1.0. The jets are required to lie 
within the region —1.0 < rf^ t < 2.2 to ensure good containment within the LAr calorimeter. 

7 For events fully contained in the detector, the total E — p z is sharply peaked at 2E e = 55 GeV. 
8 When calculating all hadronic final state quantities, particle masses are neglected. 
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Figure 3: Observed distributions of the average transverse energy flow per event around the 
jet axes in the diffractive dijet sample. A77* and A0* are the distances from the jet axes in 
pseudorapidity and azimuthal angle in the hadronic centre-of-mass frame. The jet profiles in 
77 and 4> are integrated over ±1 unit in <p and 77 respectively, (a) and (c) show the distributions 
for the backward jet in the laboratory frame, whereas (b) and (d) show those for the forward 
jet. The distributions for the simulated sample of RAPGAP events are compared with the data. 
Here, the contributions from direct photons only (dotted histograms) and from the sum of direct 
and resolved photon contributions (solid histograms) are shown. 

Transverse energies and momenta are calculated with respect to the 7*p axis. Events with either 
at least two or exactly three jets with transverse momentum pr,jet > 4 GeV are selected for the 
dijet and 3-jet samples respectively. The average resolution in Prjet ls 14%. No requirements 
are made on the presence or absence of hadronic activity beyond the jets. The final event 
selection yields 2506 dijet and 130 3-jet events. 

Fig. |3| shows the transverse energy flow around the jet axes for the dijet sample. For the 
jet profiles in 77 and <p, only transverse energies within one unit in azimuth and pseudorapidity 
are included in the plots respectively. The jet profiles for backward and forward jets are shown 
separately in Figs. |3]a,c and b,d. The data exhibit a clear back-to-back dijet structure in azimuth. 
The energy flow is well described by the RAPGAP simulation that is used to correct the data 
(solid lines). 
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Kinematic Range of 
Hadron Level Cross Sections 

4 < Q 2 < 80 GeV 2 
0.1 < y < 0.7 

xjp < 0.05 
My < 1.6 GeV 
|*j < 1.0 GeV 2 
Nj Cts > 2 or N iets = 3 
^ Jet > 4 GeV 
-3 < V* et < 



Table 1 : The kinematic range to which the dijet and 3-jet hadron level cross sections are cor- 
rected. The details of the jet finding algorithm can be found in section O. 



4.3 Cross Section Measurement 

The data are first corrected for losses at the trigger level, which occur due to the track require- 
ment. For the selected events, the trigger efficiency varies between 80 and 90%, depending 
on the kinematics. Corrections for detector acceptances and migrations between measurement 
intervals are evaluated by applying a bin-to-bin correction method using the RAPGAP program 



(see section [33]). The simulation gives an acceptable description of all relevant kinematic dis- 
tributions of the selected dijet and 3-jet events. Smearing in xjp is taken into account up to 
xjp = 0.2 by the simulation of colour-singlet exchange in RAPGAP. Migrations from xp > 0.2 
or from large values of M Y > 5 GeV are covered by a RAPGAP simulation of inclusive DIS. 
This contribution is at the level of 5% averaged over all measured bins and is concentrated at 
large xjp. An additional correction of— 6.5% ± 6.5% is applied to account for the net smearing 
about the M Y = 1.6 GeV boundary. Since only elastically scattered protons have been sim- 
ulated in RAPGAP, this correction is evaluated using the proton dissociation simulation in the 
DIFFVM [|53|] Monte Carlo model. A further correction of +5.5% ± 1.4% takes into account 
diffractive events rejected due to fluctuations in the noise level in the FMD. This correction is 
determined using randomly triggered events. QED radiative corrections are of the order of 5%. 
The bin purities and stabilities^ are typically of the order of 50 to 60% and it is ensured that they 
exceed 30% for every measured data point. 

The corrected cross sections are defined in a model independent manner, whereby the sys- 
tems X and Y are separated by the largest gap in rapidity among the hadrons in the centre- 
of-mass frame (fig. |l[). The ep cross sections are corrected to the hadron level and are quoted 
at the Born level. The kinematic range in which the cross sections are measured is fully speci- 
fied in tab. U The measured range of jet pseudorapidities in the hadronic centre-of-mass frame 
—3 < r]* et < approximately matches the —1 < rf-ft < 2.2 cut for the selected events. No rj mSLX 
or similar cuts are imposed in the definition of the measured cross sections. 



9 'Bin purity' is defined as the fraction of simulated events reconstructed in a bin that are also generated in that 
bin. 'Stability' is the fraction of events generated in a bin that are also reconstructed in that bin. 



11 



4.4 Analysis of Systematic Uncertainties 

The following sources of uncertainty contribute to the systematic errors on the measured cross 



sections. The uncertainties associated with detector understanding (see section |4.1| ) are as fol- 
lows. 



1 . The uncertainties in the hadronic calibrations of the LAr and SPACAL calorimeters mainly 
influence the measured values of firjet an ^ Mx- The resulting uncertainties in the cross 
sections are up to 10% (with a mean value of 5%) for the LAr and 0.5% for the SPACAL. 

2. The uncertainties in the E' e and 6' e measurements propagate into the reconstruction of 
Q 2 , y and W and the definition of the ^*p axis for the boost into the ^*p centre-of-mass 
frame. The uncertainty in 9' e leads to a systematic error of 1% to 2%. The uncertainty in 
E' e results in a systematic error between 1% and 5%, depending on the kinematics. 

3. The uncertainty in the fraction of energy of the reconstructed hadronic objects carried by 
tracks is 3%, leading to a systematic error in the range 1% to 5%. 

4. The uncertainties in the determinations of the trigger efficiency and the ep luminosity 
affect the total normalisation by 5% and 2% respectively. 

5. There is an uncertainty of 25% in the fraction of events lost due to noise in the FMD, 
which translates into a 1.4% normalisation error on the measured cross sections. 



The Monte Carlo modelling of the data gives rise to the following uncertainties. 



6. The uncertainty in the number of events migrating into the sample from x F > 0.2 or 
My > 5 GeV is estimated as 25%, leading to a systematic error between 1% and 3%, 
with the biggest values at large xp. 

7. A 6.5% uncertainty arises from the correction for smearing about the My limit of the 
measurement. It is estimated by variations of: (a) the ratio of elastic proton to proton 
dissociation cross sections in DIFFVM between 1:2 and 2:1; (b) the generated My dis- 
tribution within 1/My 0±0 " 3 ; (c) the t dependencies in the proton dissociation simulation 
by changing the slope parameter by ±1 GeV~ 2 and (d) the simulated efficiencies of the 
forward detectors FMD and PRT by ±4% and ±25% respectively. 

8. The uncertainty arising from the QED radiative corrections is typically 5%, originating 
from the limited statistics of the Monte Carlo event samples. 

9. The use of different approximations for higher order QCD diagrams (the parton shower 
(MEPS) model or the colour dipole (CDM) approach) leads to a 3% uncertainty in the 
cross sections. 
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10. The model dependence of the acceptance and migration corrections obtained from the 
RAPGAP simulation is estimated by varying the shapes of kinematic distributions in the 
simulations beyond the limits imposed by previous measurements or the present data. 
This is done by reweighting (a) the zjp distribution by z^' 2 and (1 — zjp) ±0 - 2 ; (b) the px 
distribution by (l/p T ) ±0 ' 5 ; (c) the xp distribution by (l/xjp) ±0 2 ; (d) the t distribution 
by e ±2t and (e) the rf-f t distribution to that observed in the data. The resulting systematic 
uncertainties range between 6% and 13%, the largest contributions originating from (c) 
and (e). 

1 1 . The lower p^-cut-off chosen to avoid collinear divergences in the leading order QCD 
matrix elements in RAPGAP is relatively high (p T > 9 GeV 2 ) with respect to the experi- 
mental cut of PTjet > 16 GeV 2 . Studying the dependence on the cut-off value results in 
an additional uncertainty of 5%. 



Most of the systematic uncertainties are not strongly correlated between data points. The 
total systematic error has been evaluated for each data point by adding all individual systematic 
errors in quadrature. It ranges between 15 and 30% and for most data points is significantly 
larger than the statistical uncertainty. 



5 Results 



In this section, the measured hadron level differential cross sections are presented for the kine- 
matic region specified in tab. [j]. The cross sections are shown graphically in Figs. |5Hl2|. In all 
figures, the inner error bars correspond to the statistical error, whilst the outer error bars rep- 
resent the quadratic sum of the statistical and systematic errors. The numerical values of the 
measured cross sections can be found in Tabs. 0-[7| The quoted differential cross sections are 
average values over the intervals specified in the tables. 



5.1 General Properties of the Dijet Data 

In this section, general features of the data are discussed, referring to Figs. [f|-[7| The model 
predictions^] which are also shown in these figures are discussed in sections |5]2] and |53| . 

In fig. ||a, the uncorrected average transverse energy flow per event for the dijet sample is 
shown as a function of the pseudorapidity in the rest frame of the X systemQ. Positive values 
of correspond to the pomeron hemisphere, negative values to the photon hemisphere. Both 
the total energy flow and the energy flow from particles outside the two leading jets are shown. 
The data exhibit considerable hadronic energy not associated with the jets. This additional 
energy is distributed in both hemispheres with some preference for the pomeron hemisphere. In 

Software to produce predictions for the measured cross sections using any hadron level ep Monte Carlo model 
is available in the HZTOOL package [^6|] . 

"This frame can be interpreted as the j*1P centre-of-mass frame. In this context, 'IP' or 'pomeron' is used 
synonymously with 'colourless exchange'. 
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Figure 4: (a) The uncorrected distribution of the average transverse energy per event for the 
diffractive dijet sample as a function of the pseudorapidity in the centre-of-mass frame of the 
X system. Distributions are shown both for all final state particles (solid points) and for only 
those particles which do not belong to the two highest pr jets (open points). The prediction 
of the RAPGAP simulations for direct and for direct plus resolved virtual photon contributions 
are also shown, (b) The uncorrected correlation between the squared invariant mass of the X 
system M\ and the squared dijet invariant mass M 2 2 for the diffractive dijet sample. The dotted 
line corresponds to M\ = Mf 2 . 



order to examine the sharing of energy within the X system on an event-by-event basis, fig. ||b 
shows the uncorrected correlation between the squared dijet invariant mass M 2 2 and the squared 
total diffractive mass M\ []25|]. M i2 is calculated from the massless jet 4-vectors. Except for a 



small subset of the events at low Mx, only a fraction of the available energy of the X system is 
contained in the dijet system, such that M\i is considerably smaller than Mx on average. 

Figs. ^] and |5| present differential dijet cross sections as functions of the following observ- 
ables: the photon virtuality Q 2 ; the mean dijet transverse momentum PT.j e ts> defined as 

Pt, jets 2 

the Yp invariant mass W; the mean dijet pseudorapidity in the laboratory frame (r])^ ets , defined 
as 

and the logarithms of the xp and (3 variables. The Q 2 and p* T ^ ets distributions are steeply 
falling. Due to the selection of events with Q 2 > 4 GeV 2 andp^ 2 ets > 16 GeV 2 , the relation 
Prjets > Q 2 holds for the bulk of the data. As can be seen in fig. ^p, the W range of the selected 
events is approximately 90 < W < 260 GeV. The xp distribution shows a rising behaviour 
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Figure 5: Diffractive dijet cross sections as a function of (a) the photon virtuality Q 2 , (b) the 
mean transverse jet momentum jets m tne 1*P centre-of-mass frame, (c) the 7*p invariant 



mass W and (d) the mean jet pseudorapidity in the laboratory frame (f]) l ^ ts . Also shown are the 
predictions from a resolved (partonic) pomeron model with gluon dominated pomeron parton 
distributions as obtained from the QCD analysis of by HI [TJ]. The results, using both the 
'fit 2' ('fiat gluon') and 'fit 3' ('peaked gluon') parton distributions for the pomeron, are shown 
evolved to a scale /j 2 = Q 2 + p\. Resolved virtual photon contributions are added according 
to the SaS-2D parameterisation [|48|]. The prediction based on 'fit 2' is also shown where the 
colour dipole approach (CDM) for higher order QCD effects is used in place of parton showers 
(MEPS). 
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Figure 6: Differential diffractive dijet cross sections as functions of (a) \ogxp and (b) log (3. 
The solid curves represent the predictions of the resolved pomeron model ('fit 2') as described 
in the text with direct and resolved photon contributions. For the log xp distribution, the con- 
tribution from sub-leading reggeon exchange is indicated by the hatched area. The dashed and 
dashed-dotted histograms correspond to the cross section predictions where the value of the 
pomeron intercept ap(0) in the model was changed from the default value of 1.20 to 1.08 and 
1.40 respectively. For this figure, all model predictions have been scaled to the integrated cross 
section in the data. For the log (3 distribution, the prediction using the 'fit 3' parton distributions 
is also shown and the range covered by the inclusive HI measurement of F® is indicated. 



from the lowest accessible values of ~ 0.003 up to the cut value of 0.05. For kinematic reasons, 
the dijet measurement is dominated by larger x p values than is the case for inclusive diffractive 
measurements. The j3 range covered by the measurement extends down to almost 10~ 3 , lower 
than accessed so far in measurements of ■ The shapes of the measured cross sections are 
generally well described by the RAPGAP simulation used to correct the data (solid histograms), 
except for the (r]) l ^f ts distribution, which shows that on average the measured jets have slightly 
larger pseudorapidities than is predicted by the simulations. 

In fig. [7], the cross section is shown differentially in z F ets \ which is calculated from 

(jets) _ Q 2 + M±2 n ~ 

Z]P ~ Q* + M% ■ ( } 

Monte Carlo studies show that the resolution in z F ets ^ is approximately 25% and that there is a 
good correlation between z P ets ^ and the value of z$> as defined at the parton level in eq. [5|. In 
loose terms, the z F ets ^ observable measures the fraction of the hadronic final state energy of the 
X system which is contained in the two jets. The measured z%^ ts ^ distribution is largest around 
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Figure 7: The diffractive dijet cross section as a function of z^ ts ^ . The same data are compared 
to predictions of resolved pomeron models, where either (a) fi 2 = Q 2 + p\ or (b) fj 2 = p\ are 
used as renormalisation and factorisation scales. In (a), the 'fit 2' (or 'fiat gluon') and 'fit 3' 
(or 'peaked gluon') parameterisations based on the HI leading order QCD fits to [|]] are 

shown. Direct and resolved 7* contributions are both included. The size of the resolved 7* 
contribution in 'fit 2' is indicated by the shaded histogram. In (b), where only the direct 7* 



contributions are shown, the preferred solution 'ACTW fit D' of the fits from Q58D is shown in 
addition to the HI fits. The corresponding gluon distributions, evolved to the mean value of the 
respective scale used and normalised such that the pomeron flux fjp/ p (x]p = 0.003, t = 0) is 
unity, are shown above the predictions. 



0.2 and thus confirms the observation from fig. |] that the total energy of the X system is typi- 
cally much larger than that contained in the jets. Diffractively scattered qq photon fluctuations 
(see section |22|) satisfy z jp = 1 at the parton level, but can be smeared to z^ 1 ^ values as low as 
0.6 because of fragmentation and jet resolution effects. Even taking this smearing into account, 
the Zp £ts ^ distribution implies the dominance of qqg over qq scattering in the proton rest frame 
picture. 



5.2 Interpretation within a Partonic Pomeron Model 

In this section, the data are compared with the Ingelman-Schlein model (section |3~2[), using the 



RAPGAP Monte Carlo model with various sets of pomeron parton distributions. In all cases un- 
less otherwise stated, the RAPGAP predictions shown use the parton shower approximation to 
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higher order diagrams (MEPS) and a contribution from resolved virtual photons is included, as 
described in section j3~3j It has been shown in an HI measurement of inclusive dijet production 
for similar ranges in Q 2 and Vt jets [0 mat including resolved photon contributions improves 
the description of the data by leading order Monte Carlo Models in the region Pr%ts > Q 2 - It 
is thus reasonable to expect a similar contribution in diffraction. 



5.2.1 Diffractive Gluon Distribution 

Pomeron parton densities dominated by gluons have proved successful in describing not only 
inclusive measurements of the diffractive structure function [0-0], but also more exclusive 
hadronic final state analyses [Q-|Sp. By contrast, pomeron parton distributions dominated by 
quarks (e.g. 'fit 1' from [Q]) do not describe the data [[3],j5|,|8|]. In particular, they lead to signif- 
icantly smaller predicted dijet electroproduction cross sections than were obtained in previous 
measurements [§]. The free parameters of the Ingelman-Schlein model to which dijet produc- 
tion is most sensitive are the pomeron gluon distribution gp(z, /i 2 ) and the pomeron intercept 
ajp(O). The sub-leading reggeon contribution and the pomeron quark distribution are better 
constrained by inclusive colour singlet exchange measurements ||3Ll57|]. 

Predictions based on two sets of pomeron parton distributions obtained from the leading or- 
der DGLAP analysis of F° {3) from HI in []3p are compared with the data in Figs. |5|, The 'flat 
gluon' or 'fit 2' parameterisation gives a very good description of all differential distributions, 
except for da / d{r]) l ^ ts . The predictions based on the 'peaked gluon' or 'fit 3' parameterisation 
in Figs. 0, |^ are also in fair agreement with the data, though the description is somewhat poorer 
than that from 'fit 2'. If the colour dipole approximation (CDM) to higher order QCD effects is 
used instead of parton showers (MEPS), the predicted dijet cross sections increase in normali- 
sation by approximately 15% (fig. |5|). The shapes of the predicted distributions, including that 
of z% ets \ are not significantly affected. 

The cross section differential in z% (fig. 0) is also compared with predictions from dif- 
ferent sets of pomeron parton distributions. Fig. [7|a shows the predictions based on the par- 
tons extracted in 'fit 2' and 'fit 3' of 0]. The parton distributions are evaluated at a scaleQ 
/i 2 = Q 2 + pip. The contribution of quark induced processes in the predictions is small. The 
fraction of the cross section ascribed to resolved virtual photons, which is shown separately for 
'fit 2' in fig. [7|a, is also small and is concentrated at low z% ets \ The same is true for the reggeon 
contribution (not shown). The predictions based on the 'fiat gluon' or 'fit 2' parton densities are 
in very good agreement with the data. The 'peaked gluon' or 'fit 3' parameterisation leads to 
an overestimate of the dijet cross section at high values of z P ets \ The gluon distributions from 
which the predictions are derived are shown above the data at jj 2 = 42 GeV 2 , representing the 
mean value of Q 2 + Prfjets f° r me selected events. The difference in shape between the gluon 
distributions and the hadron level predictions reflects the kinematic range of the measurement 
(tab. [[]). The dijet data are highly sensitive to the shape of the gluon distribution, which is poorly 
constrained by the measurements. This is especially the case in the region of large mo- 

mentum fractions (zjp or (3), since data with j3 > 0.65 were excluded from the DGLAP analysis 
ofF 2 D{3 \ 

12 Alternative reasonable choices of scale such as Q 2 + Ap^ make only small differences to the Monte Carlo 
predictions. 
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In fig. ^]b, the same data are compared with the models where p\ was chosen as the renor- 
malisation and factorisation scale and only direct photon contributions are included. The level 
of agreement between the data and the simulations based on the HI fits is similar to that in 
fig. 0a. Also shown is a prediction based on the best combined fit in [ |58"1 ] to HI and ZEUS 



F2 data and ZEUS diffractive dijet photoproduction dataF} Due to the different shape and 
normalisation of the gluon distribution in this parameterisation, the agreement with the dijet 
data is significantly poorer than is the case for the two HI fits. 

In general, the close agreement between the 'fit 2' and 'fit 3' parameterisations and the 
data can be interpreted as support for factorisable pomeron parton distributions in DIS, strongly 
dominated by gluons with a momentum distribution relatively flat in zp. 



5.2.2 Scale Dependence, Regge Factorisation and Pomeron Intercept 

In the following, some basic assumptions of the resolved pomeron model are tested, namely the 
evolution of the parton distributions with scale, Regge factorisation and the universality of the 
pomeron intercept. 

Fig. |]a shows the cross section differential in z%^ ts ^ in four intervals of the scale ji 2 = 
Q 2 + Pt- Even in this double differential view, the 'fit 2' resolved pomeron model with parton 
densities evolving according to the DGLAP equations gives a very good description of the data. 
The 'peaked gluon' solution overestimates the cross section at high z P ets ^ in all regions of /x 2 . 

In fig. |8|b, the data are used to test Regge factorisation (eq. [7]). The cross section differential 
in z% ets ^ is measured in four intervals of xp. A substantial dependence of the shape of the z^ 1 ^ 
distribution on xp is observed. This is dominantly a kinematic effect, since xjp and z% ets ^ are 
connected via the relation xp ■ z P ets ^ = x^ ets \ where x^ ets ^ is the proton momentum fraction 
which enters the hard process. The range in x^ ets ^ is approximately fixed by the kinematic range 
of the measurement. Again, the factorising resolved pomeron model describes the distributions 
well. Thus, at the present level of precision, the data are compatible with Regge factorisation. 
There is little freedom to change the pomeron intercept «jp(0) and compensate this by adjusting 
the gluon distribution. Fast variations of ap(0) with zp are also incompatible with the data. 

The value of ajp(O) controls the energy or x F dependence of the cross section. In the 
predictions of the resolved pomeron model shown in Figs. |5|-|8], a value of ap(0) = 1.2 is 
used, as obtained in the HI analysis of [fSp. Since this value of ajp(O) is larger than 

that describing soft interactions, it is interesting to investigate whether further variation takes 
place with the additional hard scale introduced in the dijet sample. In fig. the effect on the 
shape of the predicted cross section differential in x F is investigated when ap(0) is varied. As 
examples, the predictions with ap(0) = 1.08 ('soft pomeron') and ap(0) = 1.4 (approximate 
leading order 'BFKL pomeron' [|3~51]) are shown. All predictions have been scaled to the total 
cross section in the data. The xp dependence of the data requires a value for ap(0) close to 
1.2. The values of 1.08 and 1.4 result in xp dependences which are steeper or flatter than the 

13 In this parameterisation, the pomeron intercept is set to ap(Q) — 1.19. 
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Figure 8: Diffractive dijet cross sections as a function of z% ets \ shown in four intervals of (a) 
the scale /j 2 = Q 2 + "p\ and (b) log xjp. The data are compared to the resolved pomeron model 
based on the two fits to from HI, including both direct and resolved 7* contributions. 



data respectively. Making a fit for ajp(O) to the shape of the xp cross section, assuming a flux 
of the form given in eq. |9|, yields a value of 

ajp(O) = 1.17 ± 0.03 (stat.) ± 0.06 (syst.) t° f 7 (model) . 

The model dependence uncertainty is evaluated by varying the resolved photon and the reggeon 
contributions in the model by ±50% each, changing the pomeron gluon distribution within 
the range allowed by the measured z% ets ^ distribution, varying the assumed a'jp within 0.26 ± 
0.26 GeV~ 2 and varying bp between 2 GeV~ 2 and 8 GeV -2 . The effects of NLO corrections 
and possible pomeron-reggeon interference have not been studied. The extracted value of «jp(0) 
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is compatible with that obtained from inclusive diffraction in a similar Q 2 region, despite the 
fact that the jets introduce an additional hard scale. 



5.3 Energy Flow in the Photon Hemisphere and Resolved Virtual Photons 

As can be seen from Figs. |5|-|8], the data are well described by the resolved pomeron model, 
where a contribution from resolved virtual photons is included as described in section 1331. In 



this section, two observables are studied which are particularly suited to the interpretation of 
the data in terms of direct and resolved photon contributions. 

As in the case of real photoproduction analyses (see e.g. []59[j), a quantity x 7 is defined as 
the fraction of the photon momentum which enters the hard scattering. If the 4-vector of the 
parton from the photon entering the hard scattering is labelled u, then 



P-u 
P-q 



Xy = — . (16) 



Direct photon events satisfy x 7 = 1 by definition. Events where the photon is resolved have 
Xj < 1. At the hadron level, an observable x^ ets ^ can be constructed by measuring the ratio of 
the summed E — p z of the two jets to the total E — p z : 



x Uets) = ^etsE Pz 
Hall E-P 



The observable x^ ets "> correlates well with the parton level x 7 and is reconstructed with a resolu- 
tion of approximately 12% relative to the hadron level definition. The cross section differential 
in x^ ets ^ is shown in fig. [5]a. The distribution is peaked at values around 1 but there is also a 
sizeable cross section at lower x^ ets ' values. The prediction of the resolved pomeron model 
with only direct photon contributions describes the high x^ ets ^ region, but lies significantly be- 
low the data at low values of x^ ets \ The prediction is non-zero in this region only because of 
migrations from the parton level value of x 7 to the hadron level quantity x^ ets \ If the contri- 
bution from resolved photons is included, a much improved description of the data is achieved. 
The total predicted dijet cross section then increases by 17%. 

The part of the hadronic final state not associated to the two highest p* T jets is best studied in 



the 7*iP centre-of-mass frame (see section pT|). Hadronic final state particle production outside 
the two highest p* T jets can originate from jet resolution effects, possible photon and pomeron 
remnants or from higher order QCD diagrams. In order to further investigate the energy in the 
photon hemisphere, a new observable E^p m is constructed. E$) m is defined as the energy sum 
of all final state hadrons in the photon hemisphere (rf < 0) which lie outside the two highest 
p* T jet cones. The cross section is shown differentially in EyJ m in fig. |9|b. The distribution falls 
quickly as E$) m increases, indicating the dominance of direct photon scattering. The description 
at high E^lm values (corresponding to x 7 < 1) is again much improved by adding the resolved 
7* contribution. 

The presence of resolved virtual photon contributions is also suggested by the energy flow 
backward of the jets (corresponding to the photon direction) in the jet profiles (fig. Q). Similarly, 



21 



H1 Diffractive Dijets 



~ 10 3 

in ' 

X 
73 



73 



10 



10 



H1 Data (a) 

res. IP (dir.+res.y*) 
res. IP (dir.y*) 



0.2 0.4 0.6 



0.8 1 

v (jets) 



> 
0) 

O 

Q. 

2 | 

HI 

73 

73 



10 



(b) 







10 



15 

(Y) 



20 



E ™ [GeV] 

rem 



Figure 9: Differential diffractive dijet cross sections as a function of (a) x^ ets \ an estimator for 
the photon momentum fraction entering the hard scattering process, and (b) E^ m , the summed 
hadronic final state energy not belonging to the two highest p? jets in the photon hemisphere 
of the 7*^ centre-of-mass frame. The data are compared to the resolved pomeron model ('fit 
2') with and without an additional contribution from resolved virtual photons, parameterised 
according to the SaS-2D photon parton distributions. 



the transverse energy not associated with the jets in the rf < hemisphere of the 7*iP system 
(fig. f|a), is best described when the resolved photon contribution is added. Good descriptions 
of these distributions cannot be achieved by adjusting the diffractive gluon distribution. The 
resolved virtual photon contributions can be viewed as an approximation to NLO QCD diagrams 
and/or contributions without strong k T ordering. The possible presence of such effects will be 



investigated further in section 5.5 



5.4 Soft Colour Neutralisation Models 



The Soft Colour Interactions (SCI) and semiclassical models (section pA\ ) both give a reason- 
ably good description of inclusive diffraction at HERA with a small number of free parameters. 
In fig. [10|, the predictions of these models are compared with the dijet cross sections as func- 
tions of Pt jets* Mx, logxjp and z^ ts \ With the exception of the cross section differential in 
Mx, the data shown are identical to those in earlier figures. The original version of SCI gives 
a reasonable description of the shapes of the differential distributions of the dijet data, but the 
overall cross section is too low by a factor of about 2. The refined version of the SCI model, 
based on a generalised area law for string rearrangements, gives an improved description of 
F®^ 1 at low Q 2 . It also reproduces the normalisation of the dijet cross sections much better 
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Figure 10: Differential dijet cross sections as functions of (a) - eta , (b) M x , (c) logxp and 

(d) Zp Cts \ The data are compared to the original version of the Soft Colour Interaction (SCI) 
model, labelled 'SCI (original)', the prediction of the refined SCI version based on a generalised 
area law for string reconnections, labelled 'SCI (area law)', and to the semiclassical model. 



than the original version. However, the shapes of the differential distributions are not described, 
with the exception of Prjets- 

The semiclassical model gives a good description of the shapes of the distributions, but the 
total predicted dijet cross section is only around half that measured. The free parameters of the 
semiclassical model were determined using only data in the region xp < 0.01. Even at 

low xp, the predictions lie significantly below the dijet data (fig. |T0|c). It is possible that the 
inclusion of NLO terms would improve the description of the data by the semiclassical model. 
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5.5 Colour Dipole and 2-Gluon Exchange Models 



In this section, the saturation and BJLW models (section p.3| ), based on the ideas of dipole cross 
sections and 2-gluon exchange, are compared with the dijet data. Because of the nature of the 
2-gluon models, only final state parton showers are included in the simulations. A restricted 
data sample with the additional cut 

xjp < 0.01 (18) 

is studied, because the calculations were carried out under the assumption of low xjp to avoid 
contributions from secondary reggeon exchanges and ensure that the proton parton distributions 
are gluon dominated. Applying this additional restriction reduces the number of events in the 
data sample by a factor of approximately 4. 

The resolved pomeron model implies the presence of a soft pomeron remnant. The same 
is true for qqg production within the saturation model where the gluon behaves in a 'remnant- 
like' manner, due to the A; T -ordering condition imposed in the calculations. By contrast, the 
qqg calculation within the BJLW model imposes high transverse momenta on all three partons 
and is not restricted to /cy-ordered configurations. Any 'remnant' system beyond the dijets in 
this model is thus expected to have relatively large p T . To gain more insight into the properties 
of the part of the hadronic final state not belonging to the jets, a new observable Vr'rem x% 



introduced. By analogy with the definition of E]fp m (section |53J), this variable measures the 



transverse momentum of all hadronic final state particles in the pomeron hemisphere of the 
7*iP centre-of-mass frame {rf > 0) not belonging to the two highest p* T jets. 

Dijet cross sections for the region x jp < 0.01 differential in Q 2 , p* T ^ ets , z% ets ^ and Vr > )em ^ 
shown in fig. |TT} They are compared with the predictions of the saturation, BJLW and resolved 
pomeron ('fit 2') models. The saturation model is able to reproduce the shapes of the measured 
cross sections, though the overall predicted dijet rate is too low by a factor of approximately 2. 
The normalisation of the saturation model is fixed from the fit to inclusive F 2 data and by the 
assumed e 6t dependence for diffractive processes. The total predicted dijet cross section would 
increase whilst preserving a good description of F®^ 1 if the t dependence were found to be 
harder for dijet production than for inclusive diffraction. 

In the BJLW model, the contribution from qq states alone is negligibly small even at large 
values of zjp. This is in accordance with the expectation for high p T , high M x diffractive final 
states. The predicted qqg contribution is much larger. The normalisation of the BJLW model 
for qqg production can be controlled by tuning the lower cut-off on the transverse momentum 
of the gluon p™ g in the calculations. If this cut-off is set to 1.5 GeV, the total cross section 
for dijet production with xp < 0.01 is approximately correct in the model. Lowering p^ g to 
1.0 GeV leads to a prediction significantly above the measured cross section. The description 
of the shapes of the differential cross sections is reasonable apart from small discrepancies in 
the z P ^ ts ^ distribution. The differences between the predictions of the saturation and BJLW 
models may originate from the different parameterisations of T[x, h?), the different treatments 
of non-fcr-ordered configurations or from the assumed t dependences. 

The resolved pomeron model, in which the non-fc^ -ordered resolved photon contributions 

are small in the low xjp region, continues to give the best description of all observables, includ- 

(ip) (jp) 
ing the p T rem distribution. The good description of the p y T r ' em distribution by both the resolved 
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Figure 11: Diffractive dijet cross sections in the restricted kinematic range x F < 0.01, shown 
as functions of (a) Q 2 , (b) prjets^ ( c ) Zp ets ^ and (d) PT^ em , the latter denoting the summed 
transverse momentum of the final state particles not belonging to the two highest jets and 
located in the pomeron hemisphere of the 7*iP centre-of-mass frame. The data are compared 
to the saturation, BJLW and resolved pomeron ('fit 2', direct and resolved virtual photons) 
models. For the BJLW model, the contribution from qq states alone and the sum of the qq and 
qqg contributions for two different values of the p T cut-off for the gluon p^ n * are shown. 



pomeron and the BJLW models indicates that the present data are not easily able to discriminate 
between models with a soft 'remnant' and those with a third high-p T parton. 
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Figure 12: Differential cross sections for diffractive 3-jet production as functions of (a) the 3- 
jet invariant mass M123 and (b) the corresponding ^-variable z^> ^ , measuring the colourless 
exchange momentum fraction which enters the hard interaction. The data are compared with 
the resolved pomeron model with two different approaches for higher order QCD diagrams, 
the parton shower model (labelled 'MEPS') and the colour dipole approach (labelled 'CDM'). 
The 'HI fit 2' parameterisation is used and direct and resolved virtual photon contributions are 
included. The BJLW model is also shown, including qq and qqg contributions, with the cut-off 
for the gluon set to 1 .5 GeV. 



5.6 3-Jet Production 



The diffractive production of three high-p T jets as components of the X system has been inves- 
tigated. Except for the requirement on the number of jets, the analysis is identical to the dijet 
analysis, such that no requirements are made on possible hadronic activity beyond the jets. In 
fig. [12], the measured 3-jet cross sections are presented as functions of the 3-jet invariant mass 
Mi 23 and 

(3 jets) _ Q 2 + M123 n q\ 

Z]P ~ Q2 + J14 ' (iy) 

which, similarly to z^ 1 ^ for dijet events, is a measure of the fraction of the energy of the 
X system which is contained in the jets. The z% J6<s ' ) cross section is measured up to 0.8. 
With the present statistics, it is not possible to extract a cross section for the interesting region 

(3 jets') 

0.8 < z)p < 1.0, which corresponds approximately to 'exclusive' 3-jet production. The 
measured ^ ets ^ cross section demonstrates that additional hadronic activity beyond the jets is 
typically present even in the 3-jet sample. 

The data are compared with the resolved pomeron model ('fit 2'), with the hard interaction 
evaluated at a scale /x 2 = Q 2 + p\. Direct and resolved 7* contributions are included. Because 
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the leading order for 3-parton final states is 0(a 2 ), two different approximations for higher 
order QCD diagrams are considered here, the parton shower model (MEPS) and the colour 
dipole approach (CDM). The measured cross sections are well described when using CDM. 

(3 jets) 

The MEPS simulation tends to lie below the data at low Mi 23 or high z F . 

The BJLW calculation with p c ^ l g = 1.5 GeV is not able to accommodate the observed rate of 

3-jet events. The predicted cross section increases towards the high z% regime of exclusive 
3-jet production. For kinematic reasons, the 3-jet sample originates from the region x jp > 0.01, 
where contributions from the proton quark distributions and secondary exchanges, which are not 
included in the 2-gluon models, can no longer be neglected. An improvement in the predictions 
of dipole models may also come through the inclusion of higher multiplicity photon fluctuations 
such as qqgg, which have not yet been calculated. 



6 Summary and Final Remarks 

An analysis of the production of jets as components of the dissociating photon system X in the 
diffractive DIS reaction ep — > eXY has been presented for 4 < Q 2 < 80 GeV 2 , xp < 0.05, 
Pt jet > ^ GeV, My < 1.6 GeV and \t\ < 1.0 GeV 2 . The kinematic range has been extended 
to lower Q 2 and VTA&t compared to previous measurements [§] and the statistical precision is 
much improved. Cross sections for the production of three high transverse momentum jets have 
been measured for the first time in diffraction. 

The observed dijet events typically exhibit a structure where, in addition to the reconstructed 
jets, the X system contains hadronic energy with transverse momentum below the jet scale. The 
dijet invariant mass is thus generally smaller than M X - Viewed in the proton rest frame, the 
data clearly require the dominance of higher multiplicity photon fluctuations (e.g. qqg) over the 
simplest qq configuration. Considered in the proton infinite momentum frame, the data show 
that the diffractive gluon distribution is much larger than the quark distribution. 

The data can be described by a 'resolved partonic pomeron' model, with diffractive parton 
distributions extracted from F2 data. The good description from this model strongly supports 
the validity of diffractive hard scattering factorisation in DIS. The dominant contribution in the 
model arises from a diffractive exchange with factorising xp dependence ('Regge' factorisa- 
tion). A value of a P (0) = 1.17 ± 0.03 (stat.) ± 0.06 (syst.) t° f 7 (model) is obtained for the 
intercept of the leading trajectory from fits to the dijet data. The compatibility of the data with 
QCD hard scattering and Regge factorisation contrasts with the observed strong factorisation 
breaking when diffractive ep and pp data are compared [ |T3| , |5%| ]. The dijet data give the best con- 
straints to date on the pomeron gluon distribution. The data require a large fraction (80 — 90%, 
as obtained in [§]) of the pomeron momentum to be carried by gluons with a momentum dis- 
tribution which is comparatively flat in z F . Predictions derived from the 'flat gluon' (or 'fit 2') 
parameterisation in 0, with higher order QCD effects modelled using parton showers, are in 
remarkably good agreement with all aspects of the dijet data with the single exception of the 
{ r i) l je\ s dependence. The level of agreement between the resolved pomeron model and the data 
is better than that obtained from leading order predictions for inclusive ep dijet data (e.g. [jfjf]]), 
where the NLO corrections are approximately 40% in a similar region of Q 2 and Pxjets- 
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The two versions of the Soft Colour Interactions (SCI) model are not able to reproduce the 
overall dijet rate and the shapes of the differential cross sections at the same time. The similarly 
motivated semiclassical model in its present (leading order) form achieves a good description 
of the shapes of the differential distributions but underestimates the total dijet cross section. 

Models based on colour dipole cross sections and 2-gluon exchange have been compared 
with the dijet data in the restricted region xjp < 0.01. The saturation model, which takes 
only k T ordered configurations into account, describes the shapes of the jet distributions but 
underestimates the overall cross section. The normalisation of the BJLW model, in which strong 
kr ordering is not imposed, is close to the data if a cut-off for the gluon transverse momentum 
of p^ g = 1.5 GeV is chosen for the qqg contribution. The shapes of the differential distributions 
are reasonably well described. 

Strong conclusions cannot yet be drawn from the 3-jet production cross sections, because of 
the limited statistical accuracy and the kinematic restriction to large xp implied by the require- 
ment of three high p T jets. At the present level of precision, the partonic pomeron predictions 
based on the 'fit 2' parameterisation in [JJ] are in good agreement with the 3-jet cross sections, 
provided the CDM model of higher order QCD effects is used. The BJLW model is unable to 
reproduce the rate of observed 3-jet events when p^ is kept fixed at 1.5 GeV. 

In conclusion, diffractive jet production has been shown to be a powerful tool to gain in- 
sight into the underlying QCD dynamics of diffraction, in particular the role of gluons. The jet 
cross sections are sensitive to differences between phenomenological models which all give a 
reasonable description of ■ Models based on fully factorisable diffractive parton distribu- 
tions continue to be successful. Progress in calculations based on 2-gluon exchange has led to 
improved agreement with the data. 
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Dijet cross section as a function of Q'^. 


Bin 


Q 


* [GeV 


2 ] 


cr [pb/GcV^] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


4.0 




6.0 


21.4 


4.5 


16.4 


17.0 


2 


6.0 




10.0 


13.0 


4.0 


16.1 


16.6 


3 


10.0 




15.0 


6.3 


4.8 


17.0 


17.7 


4 


15.0 




20.0 


4.1 


6.2 


16.7 


17.8 


5 


20.0 




30.0 


2.3 


5.8 


16.4 


17.4 


6 


30.0 




40.0 


1.2 


8.0 


16.3 


18.2 


7 


40.0 




50.0 


0.7 


10.4 


19.7 


22.3 


8 


50.0 




60.0 


0.7 


12.5 


23.9 


27.0 


9 


60.0 




80.0 


0.4 


11.9 


29.6 


31.9 


Dijet cross section as a function of iets . 


Bin 


Pt, 


jets [GeV] 


cr [pb/GeV] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


4.0 




5.0 


74.9 


4.0 


17.8 


18.3 


2 


5.0 




6.0 


59.5 


3.3 


16.4 


16.8 


3 


6.0 




7.5 


28.8 


3.9 


17.3 


17.8 


4 


7.5 




9.0 


9.9 


7.0 


17.9 


19.2 


5 


9.0 




11.0 


3.4 


11.0 


17.7 


20.8 


6 


11.0 




14.0 


0.9 


18.9 


18.6 


26.5 


Dijet cross section as a function of (v) 1 ^- 


Bin 




a[pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


-1.00 




-0.66 


22.4 


13.5 


34.3 


36.9 


2 


-0.66 




-0.33 


68.9 


6.3 


17.7 


18.8 


3 


-0.33 




0.00 


112.8 


4.7 


15.7 


16.4 


4 


0.00 




0.33 


131.6 


4.2 


15.7 


16.3 


5 


0.33 




0.66 


127.9 


4.3 


17.5 


18.0 


6 


0.66 




1.00 


85.3 


5.1 


17.4 


18.2 




1.00 




1.50 


16.4 


6.8 


25.6 


26.5 


Dijet cross section as a function of M x - 


Bin 


M x [GeV] 


cr [pb/GeV] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


8.0 




14.0 


1.9 


10.1 


20.0 


22.4 


2 


14.0 




20.0 


7.5 


4.4 


15.1 


15.7 


3 


20.0 




30.0 


7.3 


3.2 


16.9 


17.2 


4 


30.0 




40.0 


4.5 


4.0 


17.8 


18.3 


5 


40.0 




60.0 


1.2 


6.2 


27.1 


27.8 


Dijet cross section as a function of W. 


Bin 


W [GeV] 


cr [pb/GeV] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


90.0 




115.0 


1.1 


6.6 


20.3 


21.3 


2 


115.0 




140.0 


1.4 


5.1 


18.5 


19.2 


3 


140.0 




165.0 


1.7 


4.4 


18.0 


18.5 


4 


165.0 




190.0 


1.3 


4.5 


17.7 


18.3 


5 


190.0 




215.0 


1.1 


4.7 


17.7 


18.3 


6 


215.0 




240.0 


0.9 


5.4 


17.0 


17.8 


7 


240.0 




260.0 


0.5 


10.3 


28.5 


30.3 



Table 2: Differential hadron level dijet cross sections. Here and elsewhere, the quoted differ- 
ential cross sections are average values over the specified intervals. 
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Dijet cross section as a function of log 10 x P . 


Bin 


logio x P 


a[pbj 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 
2 
3 
4 
5 
6 


-2.5 
-2.3 
-2.1 
-1.9 
-1.7 
-1.5 


-2.3 
-2.1 
-1.9 
-1.7 
-1.5 
-1.3 


7.3 
35.4 
88.2 
171.2 
269.3 
440.7 


21.8 
10.8 
6.8 
4.7 
3.6 
3.2 


28.8 
25.1 
17.5 
16.3 
16.3 
18.6 


36.1 
27.4 
18.8 
17.0 
16.7 
18.8 


Dijet cross section as a function of log 10 /3. 


Bin 


logio P 


u [pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 
2 
3 
4 

5 
6 
7 
8 


-2.8 
-2.5 
-2.2 
-1.9 
-1.6 
-1.3 
-1.1 
-0.8 


-2.5 
-2.2 
-1.9 
-1.6 
-1.3 
-1.1 
-0.8 
-0.5 


24.9 
88.3 
129.9 
152.7 
145.9 
85.0 
53.4 
13.5 


11.3 
5.6 
4.3 
3.9 
4.3 
7.0 
7.8 

17.7 


26.4 
18.1 
16.7 
17.4 
16.8 
17.5 
17.4 
29.8 


28.7 
19.0 
17.2 
17.9 
17.3 
18.8 
19.0 
34.6 


Dijet cross section as a function of z^ ets) . 


Bin 




<r[pb] 


stat. err. [%} 


syst. err. [%] 


tot. err. [%] 


1 
2 
3 
4 
5 
6 
7 
8 


0.000 - 
0.125 - 
0.250 - 
0.375 - 
0.500 - 
0.625 - 
0.750 - 
0.875 - 


0.125 
0.250 
0.375 
0.500 
0.625 
0.750 
0.875 
1.000 


269.4 
493.9 
331.3 
233.2 
174.2 
94.0 
39.8 
30.0 


5.8 
3.8 
4.2 
4.9 
5.9 
8.1 
11.7 
16.7 


23.7 
18.4 
18.6 
18.5 
16.1 
16.3 
16.3 
24.5 


24.4 
18.8 
19.1 
19.2 
17.2 
18.2 
20.0 
29.7 


Dijet cross section as a function of x^ ets) . 


Bin 




<r[pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 
2 
3 
4 
5 


0.0 
0.2 
0.4 
0.6 
0.8 


0.2 
0.4 
0.6 
0.8 
1.0 


25.4 
104.8 
153.8 
331.5 
428.3 


14.3 
6.5 
5.0 
3.6 
3.1 


35.1 
17.7 
18.1 
18.0 
16.7 


37.9 
18.9 
18.8 
18.3 
17.0 


Dijet cross section as a function of E^lln- 


Bin 


E&L [GeV] 


a [pb/GeV] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 
2 
3 
4 


0.0 
4.0 
8.0 
12.0 


4.0 
8.0 
12.0 
20.0 


38.1 
9.0 
4.0 
2.0 


2.5 
4.7 
6.7 
8.0 


17.2 
17.6 
25.4 
38.6 


17.3 
18.2 
26.3 
39.4 



Table 3: Differential hadron level dijet cross sections (continued). 
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Dijet cross section as a function of Zjp ets) for -1.5 < log 


io x p < -1-3. 




Bin 


Z IP 


CT [pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.00 




0.15 


232.9 


6.0 


29.7 


30.3 


2 


0.15 




0.30 


209.4 


5.3 


24.2 


24.8 


3 


0.30 




0.50 


85.4 


6.4 


20.8 


21.8 


4 


0.50 




0.70 


30.9 


10.4 


18.8 


21.5 


5 


0.70 




1.00 


3.4 


28.9 


47.0 


55.1 


Dijet cross section as a function of z^ ets) for -1.75 < log 10 x P < -1.5. 


Bin 




<r[pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.0 




0.2 


97.1 


6.4 


20.8 


21.8 


2 


0.2 




0.4 


134.3 


5.3 


19.0 


19.7 


3 


0.4 




0.6 


63.4 


7.1 


16.3 


17.7 


4 


0.6 




0.8 


21.8 


12.6 


16.6 


20.8 


5 


0.8 




1.0 


8.5 


25.8 


34.4 


43.0 


Dijet cross section as a function of z^ ets) for -2.0 < log 


io %p < -1.75. 




Bin 




<r[pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.00 




0.30 


37.8 


8.3 


23.4 


24.9 


2 


0.30 




0.45 


59.7 


9.2 


18.2 


20.4 


3 


0.45 




0.60 


49.1 


11.2 


19.4 


22.4 


4 


0.60 




0.80 


33.8 


10.9 


19.1 


22.0 


5 


0.80 




1.00 


9.0 


21.8 


24.8 


33.0 


Dijet cross section as a function of z\p ets) for log 10 x F < -2.0. 


Bin 


Zjp 


a [pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.00 




0.30 


5.2 


30.2 


85.2 


90.4 


2 


0.30 




0.45 


25.8 


16.2 


31.5 


35.4 


3 


0.45 




0.60 


28.0 


13.9 


21.3 


25.5 


4 


0.60 




0.80 


24.8 


12.8 


18.9 


22.9 


5 


0.80 




1.00 


11.9 


17.4 


22.0 


28.0 



Table 4: Differential hadron level dijet cross sections in four bins of log 10 xjp. 



Dijet cross section as a function of 


z { i et3) for 20 GeV 2 < Q 2 + p 2 


< 35 GeV 2 . 


Bin 


Zjp 


a[pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.0 - 


0.2 


150.5 




6.7 


32.4 


33.1 


2 


0.2 - 


0.4 


109.0 




7.3 


26.6 


27.6 


3 


0.4 - 


0.6 


45.2 




10.8 


28.8 


30.8 


4 


0.6 - 


0.8 


18.7 




16.2 


31.6 


35.5 


5 


0.8 - 


1.0 


5.9 




31.6 


54.4 


63.0 


Dijet cross section as a function of 


z\p ets> for 35 GeV 2 < Q 2 + p 2 


< 45 GeV 2 . 


Bin 


Zjp 


a[pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.0 - 


0.2 


89.7 




7.1 


25.2 


26.2 


2 


0.2 - 


0.4 


71.8 




6.9 


21.8 


22.9 


3 


0.4 - 


0.6 


39.3 




9.0 


26.1 


27.6 


4 


0.6 - 


0.8 


16.9 




14.4 


26.3 


30.0 


5 


0.8 - 


1.0 


4.3 




27.7 


26.2 


38.1 


Dijet cross section as a function of 


z\p ets> for 45 GeV 2 < Q 2 + p 2 


< 60 GeV 2 . 


Bin 


Zjp 


a [pb] 


stal 


. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.0 - 


0.2 


74.6 




7.7 


24.6 


25.8 


2 


0.2 - 


0.4 


78.0 




6.7 


24.5 


25.4 


3 


0.4 - 


0.6 


43.2 




8.6 


18.6 


20.5 


4 


0.6 - 


0.8 


14.7 




14.7 


20.1 


25.0 


5 


0.8 - 


1.0 


5.5 




23.6 


28.7 


37.2 


Dijet cross section as a function of 


Zjp 


s; for Q 2 


+ p 2 > 60 GeV 2 




Bin 


Zjp 


a [pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.0 - 


0.2 


58.7 




9.2 


25.2 


26.9 


2 


0.2 - 


0.4 


114.6 




5.8 


17.4 


18.4 


3 


0.4 - 


0.6 


73.4 




6.8 


15.6 


17.0 


4 


0.6 - 


0.8 


45.3 




9.4 


15.3 


18.0 


5 


0.8 - 


1.0 


14.4 




18.3 


22.3 


28.8 



Table 5: Differential hadron level dijet cross sections in four bins of Q 2 + p\. 
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Dijet cross section as a function of Q 2 for xp 


< 0.01. 






Bin 


[GeV^] 


a [pb/GcV^] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


4.0 - 10.0 


1.58 


9.7 


18.5 


20.9 


2 


10.0 - 20.0 


0.40 


13.9 


18.4 


23.1 


3 


20.0 - 40.0 


0.12 


17.7 


29.2 


34.1 


4 


40.0 - 80.0 


0.01 


44.7 


58.3 


73.5 


Dijet cross section as a function of '[f T jets for x F < 0.01. 


Bin 




a [pb/GeV] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


4.0 - 5.0 


8.3 


12.0 


23.9 


26.8 


2 


5.0 - 6.0 


4.7 


11.0 


21.3 


24.0 


3 


6.0 - 7.5 


2.0 


16.0 


19.1 


24.9 


4 


7.5 - 9.0 


0.3 


50.0 


43.2 


66.1 


Dijet cross section as a function of z\^ et3) tor x F < 0.01. 


Bin 




<r[pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.00 - 0.30 


5.2 


30.2 


85.2 


90.4 


2 


0.30 - 0.45 


25.8 


16.2 


31.5 


35.4 


3 


0.45 - 0.60 


28.0 


13.9 


21.3 


25.5 


4 


0.60 - 0.80 


24.8 


12.8 


18.9 


22.9 


5 


0.80 - 1.00 


11.9 


17.4 


22.0 


28.0 


( IP) 

Dijet cross section as a function of p^, r ; em for rr^p < 0.01. 


Bin 


P ( rZm [GeV] 


a [pb/GeV] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.0 - 0.5 


24.5 


10.3 


26.8 


28.8 


2 


0.5 - 1.0 


11.2 


12.7 


24.8 


27.9 


3 


1.0 - 3.0 


1.3 


16.0 


58.1 


60.3 



Table 6: Differential hadron level dijet cross sections in the restricted kinematical range xjp < 
0.01. 



3-jet cross section as a function of M123. 1 


Bin 


M123 [GeV] 


a [pb/GeV] 


stat. err. [%} 


syst. err. [%] 


tot. err. [%] 


1 


12.0 - 20.0 


0.48 


14.4 


33.0 


36.1 


2 


20.0 - 30.0 


0.43 


11.6 


23.0 


25.7 


3 


30.0 - 40.0 


0.06 


35.3 


40.5 


53.7 


3-jet cross section as a function of 2^ i ets \ 


Bin 


(3 jets) 

Zjp 


CT [pb] 


stat. err. [%] 


syst. err. [%] 


tot. err. [%] 


1 


0.2 - 0.4 


13.1 


16.2 


43.8 


46.7 


2 


0.4 - 0.6 


13.9 


13.9 


20.5 


24.7 


3 


0.6 - 0.8 


10.5 


18.6 


22.6 


29.2 



Table 7: Differential hadron level 3-Jet cross sections. 
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